Asymmetric Planetary Nebulae III 

ASP Conference Series, Vol. XXX, 2004 

M. Meixner, J. Kastner, B. Balick and N. Soker eds. 



MHD Disk Winds in PNe and pPNe 



-i— > 
O 



OS 



o 

in 



Adam Frank 

Department of Physics and Astronomy Department, University, of 
Rochester, Rochester, NY, 14-610, 



' afrank@pas.rochester.edu 

o ' 



Abstract. Winds from accretion disks have been proposed as the driv- 
ing source for precessing jets and extreme bipolar morphologies in Plan- 
, etary Nebulae (PNe) and proto-PNe (pPNe). Here we apply MHD disk 

wind models to PNe and pPNe by estimating separately the asymptotic 
MHD wind velocities and mass loss rates. We show that the resulting 
winds can recover the observed momentum and energy input rates for 
PNe and pPNe. 
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£S) , Introduction While considerable progress has been made in understanding the 

hydrodynamic shaping of elliptical and bilobed PNe (Balick & Frank 2002), the 
origin of extreme butterfly nebulae as well as jets in PNe continues to pose a 
number of problems for theorists. 
^5 , In particular a formidable problem for pPNe concerns the total momentum 

and energy in the outflows. A number of observational studies have shown that 
radiatively accelerated winds in pPNe cannot account for the high momentum 
and energy implied by CO profiles (Bujarrabal et aZ.2001). Thus the total flow 
momentum n is such that n >> (L*/c)Ai, where L* is the stellar luminosity 
emitted during the pPNe outflow expansion lifetime At. In light of these results 
both the launching and collimation of winds in pPNe becomes problematic. Note 
that while both the dominant hydrodynamic theory for shaping PNe, (the Gen- 
eralized Interacting Stellar Winds model: Balick 1987; Icke 1988), and models 
/\ \ invoking an initially weak toroidal magnetic field (Garcia-Segura et al. 1999) can 

■ produce jets neither can account for the momentum excesses in pPNe. 

Thus there remains considerable uncertainty about the processes which pro- 
duce collimated jets/outflows in pPNe and PNe. Other systems which produce 
jets such as YSOs, AGN and micro-quasars have been modelled via a combina- 
tion of magnetic and centrifugal forces from accretion disks (Konigl & Pudritz 
2000). The success of these Magneto- centrifugal Launching (MCL) models is 
such that it is worthwhile considering if such a scenario can be applied to PNe 
and pPNe. While models such as Morris 1987, Soker & Livio 1994 and Soker 
& Rapport 2001 have relied heavily on collimated winds from disks these works 
did not specify how such winds are launched or collimated. Thus application 
of MCL disk wind models to PNe and pPNe would close an important gap in 
building a new MHD disk wind paradigm for these systems (Blackman, Welch 
& Frank (2001): BFW01, Blackman et al. 2001). 

Here we summarize new calculations from Frank & Blackman 2004 who 
derive scaling relations from the equations for MCL and separately estimate 
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MHD disk wind mass outflow rates and asymptotic outflow velocities. These 
relations are then applied to PNe and pPNe. In these calculations a disk dynamo 
is invoked to produce the requisite fields. 

Disk Accretion Rate in PNe: In order to produce a more detailed com- 
parison of MHD disk winds with PNe it is necessary to have a model for PNe 
accretion disks. It is unlikely that an accretion disk could survive the long main 
sequence lifetime of a PN central star. Thus, unlike YSOs and AGN, accre- 
tion disks in PNe systems must form via binary interactions. Disks may form 
around secondaries via Roche lobe overflow or accretion of the dense AGB wind 
(Mastrodemos k, Morris 1998). Such systems would be similar to symbiotic 
stars. 

Mastrodemos & Morris (1998) found steady accretion disks could form 
around a white dwarf companion orbiting a AGB star with M a gb ~ 10 -5 M yr _1 . 
The ratio they found in their models of Md/M a gb ~ -05 — .005 is consistent 
with expectations from basic theory. 

Accretion disks could also form around the primary after CE evolution and 
disruption of the secondary star (Reyes- Ruiz & Lopez 1999: RRL99). This 
model implies a finite lifetime for the disk as the mass reservoir of the disrupted 
companion is slowly drained onto the primary. A description of disk formation in 
PNe has been given in RRL99 who found that systems with a primary consisting 
of an evolved AGB star with mass M* « 2.6 — 3.6 M , a low mass secondary 
(< 0.08 M ) and an initial binary separation of < 2OOi? may produce disks. 

RRL99 also found the disk accretion rate to evolve in time with in a power- 

/ f \-5/4 

law manner. Md = Mdo ( Tj/r J ^© ^ x Typical values of the scale is M^ = 
10~ 3 M yr" 1 . 

Mass Outflow Rate and Wind Speed from MCL Theory: The basic 
physics of magneto-centrifugal launching of winds and jets is well studied when 
a magnetic field distribution is imposed on the disk. In Frank Sz Blackman 2004 
it was shown how to combine Poynting flux driven outflows with asymptotic 
wind solutions and mean field dynamo theory to estimate the asymptotic wind 
speed and the outflow accretion rates. 

Magneto-centrifugal launching is a means of converting gravitational bind- 
ing energy in an accreting source into kinetic energy of an outflowing wind. The 
magnetic fields act as a drive belt to extract angular momentum from the an- 
choring rotator and launch the wind. The magnetic luminosity, or equivalently, 
the maximum magnetic power available for a wind can be obtained from the 
integrated Poynting flux (BFW01). 



where vq is the disk inner radius and Ba = B^ ~ B p at the Alfven surface 

More work is required to estimate the mass outflow rate and outflow speed 
separately. The MCL problem requires the construction of solutions for a steady, 
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ideal, isothermal magnetohydro dynamic flow. The isothermal assumption elim- 
inates the need for solving the energy equation, but more complex assumptions 
can be used, e.g. a polytropic law. Blackman & Frank 2004 began with an ex- 
pression for the magnetic Bernoulli constant which connects the wind at infinity 
to that at the footpoints of the disk 

^{u 2 p + n 2 r 2 ) + $ + w + n o {n o r\ - VLr 2 ) = U(a) = const(a). (3) 

where u p is the polodial speed at infinity, r& is the Alfven radius where u p = ua- 
This equation can then be combined with the momentum equation for the flow 
with an effective potential for a cold plasma parcel tied to a rotating field line 
(Blandford & Payne (1982)) 

u .V^- S ,^(f^-^). <4) 

Beginning with these two equations one can estimate the properties of the wind 
at large distances from the disk source if the field strength in the disk can be 
derived. 

Magnetic fields may form in these disks via dynamo processes. The topology 
of such a field, (the ratio of poloidal B p and toroidal B$ in the disk) and its 
subsequent value in the coronae remains a subject of considerable discourse. 
Frank & Blackman 2004 assumed a dynamo driven field in the disk that led 
to a primarily polodial field in a magneto-hydrostatic disk corona where the 
wind launches. From the expressions derived for the disk field along with those 
estimated from the governing MHD equations Frank & Frank 2004 provided 
estimates for the wind properties. We quote these below, 

M w ~ 0.1a ss ^M d , (5) 
ho 

■Uoo ~ 2.1O r . (6) 

L mag ~0.22a ss ^M d n 2 rl (7) 

II ~ Lmag/uoo ~ 0.Ua ss —M d tt r . (8) 

ho 

where a ss is a dimensionless parameter associated with the disk viscosity 
and ho is the disk scale height. The third and fourth expressions give estimates 
of the rate in which energy and momentum are input by the MCL disk wind 
into the ambient medium. 

Disk Winds Models for PNe and pPNe: For "classic" PNe, a total mass of 
M pn ~ .1 M.q must be accelerated to velocities of u pn ~ 40 km/ s in a timescale 
of order At pn « 10000 y. This gives II = M pn u pn / At pn ss 10 27 g cm/s 2 and 
E = M pn u 2 pn /At pn w 10 34 erg/s. 

For pPNe Bujarrabal et al. (2001) found high total outflow momentum 
10 36 < 11/ (g cm s~ l ) < 10 40 and total outflow energy 10 41 < E/(erg s' 1 ) < 
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10 47 . These values can be converted into momentum and energy injection rates 
using an assumed injection or "acceleration" timescale At, II = II/ At, L = 
E/At. 

The question which arises is: can these energy and momentum budgets be 
met with disk wind models. In what follows we use the relations from Frank & 
Blackman 2004 quoted above and assume that the inner edge of the disk extends 
to the stellar surface and use ro = r\ = r*. 

PNe Solutions: In this case we assume that the star which produces the jet is a 
proto-WD with an AGB companion (Soker h Rappaport 2001). Thus accretion 
rates of « 10~ 6 Mqj/^ 1 are reasonable. To evaluate the expressions above 
we choose M s = 0.6 M and a disk with a ss = 0.1 and vq/Kq = 10. 

For PNe central star parameters (T* = 10 5 K, L* = 5000 L such that 
ri = 1.64 x 10 10 cm) we find the following conditions for the wind from the 
equations above 

M w = Ixr'M,^ ^^ ) <») 

»~ = 1 - 25xl03to! / S (^fc) I/2 (^) _1/2 - < 10 > 

(11) 



Thus using typical conditions for PNe central stars, the scaling relations de- 
rived from the MHD equations yield disk wind parameters well matched with 
observations. 

pPNe Solutions: While the mass loss rates and velocities are known for PNe 
winds the situation for pPNe is not as clear. In general what is observed in pPNe 
is the total mass in the outflows. Velocities are also uncertain as only properties 
of swept-up material may be directly determined. 

We assume a post-AGB star with mass M s = 0.6 M , T* = 10,000 K 
and L* = 5 x 10 3 L which, assuming a blackbody, yields a radius of r* = 
1.6 x 10 12 cm = 23 R Q . Note that such a star has an escape velocity of 
u esc = 98 km s^ 1 

Achieving the high momentum input rates observed in pPNe via MCL disk 
wind models will necessitate high accretion rates. We use an accretion rate of 
Md = 1 x 10~ 4 M yr^which is the 200 year average of that found by RRL99 
for their case A Once again we choose a disk with a ss = 0.1 and ro/ho = 10. 

Assuming ro = r* along with the other parameter values given above, the 
key wind quantities M w , u w , L rn = L w , E are, 

M w = 1. x HT 5 M yr" 1 ( -] (12) 

° y \ 10~ 4 M yr^ 1 J V ; 

/ tut \ 1/2 / o \ -1/2 

»~ ^ 14 " m K^fc) ink) (13) 



.34 -1 / M d \ ( M* \ ( Ri 



L m ~ 6.7 x W' n erg s" 1 ^ T — — (14) 

y I lO- 4 M yr-iH.6M A23R ; 1 ' 
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Note that L w At ~ 10 44 , a value in the middle of the range found by Bujarrabal 
et a/.2001. Note also that the solution above has « 1.5u esc . Since Uoo u esc 
the higher velocity outflows seen in some pPNe would require more disks around 
more compact central sources. 

Given a model for the temporal history of the disk accretion, the total 
energy and momentum for the outflows can be found. Replacing with Mdif) 
from RRL99 and integrating gives 



These results show that the MCL disk wind models can achieve both energy 
and momentum injection rates as well as the total energy and momentum re- 
quired to account for many pPNe described by Bujarrabal et al (2001). The 
total energy and momenta budgets we find from these solutions fall well within 
the range of pPNe outflows with momentum excesses. Taken together with our 
previous calculations for "classic" PNe winds, these results confirm the predic- 
tions of BFW01 that magnetized disk winds can account for much of the outflow 
phenomena associated collimated outflows in the late stage of stellar evolution. 

The results above indicate that collimated flows which form from transient 
disks in the pPNe stage will appear as dense knots in mature PNe flows. This 
may also serve to explain the presence of so-called FLIERS (Fast Low Ionization 
Emission Regions) seen in some PNe. The mass loss rate in the winds derived 
above rapidly decrease with time. Thus the bulk of the jet's mass will lie near 
its head. As the material in the disk is accreted onto the star the jet will 
eventually shut-off leaving the dense knot to continue its propagation through 
the surrounding slow wind. 

When the star makes its transition to a hot central star of a PNe its fast, 
tenuous spherical wind sweeps up a shell of the slow AGB wind material. The 
shell's expansion speed will typically be of order 40 km/ s and it will not catch 
up to the head of the jet. Thus during the PNe phase the jet head will appear as 
a dense, fast moving knot which should lie outside the PNe wind blown bubble. 
We note that masses of FLIERs are estimated to be of order 10" 4 — 10 -5 M Q 
which is reasonable for the models presented above. FLIER velocities can be 
lower than the « 100 km/s calculated above but deceleration of the jet head will 
occur via interaction with the environment. We note also that hydrodynamic 
simulations of PNe jets in which the jet ram pressure decreases in time (as 
would occur for our model) show characteristic patterns of backward pointing 
bow-shocks (apex pointing back towards to the star). If such results are robust, 
the jets produced by disk winds in our scenario above may also yield similar 
morphologies. 

Discussion and Conclusions Our results for pPNe show that momentum ex- 
cesses need not occur for outflows driven by MCL winds. While this is encour- 
aging in terms of finding a mechanism for driving pPNe outflows the solutions 
require fairly high accretion rates (> 10~ 5 M Q yr _1 ). It is not clear if such 
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conditions can be achieved with the frequency required by observations. While 
solutions of RRL99 yield accretion rates and time dependencies which lead to 
the correct outflow momenta and energetics, their models place fairly stringent 
limitations on the nature of the binaries that form disks from disrupted com- 
panions. If accretion onto undetected compact orbiting companions is invoked 
Soker & Rappaport (2001) then higher values of may not be required. 

We note that a robust prediction of our models is the ratio of wind mass 
loss rate to accretion rate, i.e. M w /M a ~ .1. This is true for most MCL disk 
wind models and can be seen as a target prediction which can be explored 
observationally. 

This paper comprises a step beyond Blackman Frank Welch (2001) in estab- 
lishing the efficacy of MHD paradigms for pPNe/PNe in which strong magnetic 
fields play a role in both launching and collimating the flows. It is also worth 
noting that recent simulation results by Matt, Blackman & Frank 2004 con- 
firm that the exposed rapidly rotating magnetic core model can produce well 
collimated outflows. 
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